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ABSTRACT 


An existing code for calculating axial turbine performance 
using multiple stream surfaces was modified and made to run 
on the equivalent of an HP-1000 computer system. Calculations 
were made for the geometry of a 485 horsepower dual-discharge 
air-drive turbine for both on and off-design conditions. The 
results were compared with available data obtained at off- 
design speeds. Agreement of the flow rate and horsepower to 


within 5% was obtained. 
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I. INTRODUCTION 


A. DESCRIPTION OF THE TRANSONIC COMPRESSOR TEST RIG 

The Transonic Compressor Test Rig at the Turbopropulsion 
Laboratory (TPL) of the Naval Postgraduate School is shown 
schematically in Fig. 1 and consists of the following major 
components: 

Pawar drave turbine. 

2. Air supply system. 

3. Associated piping including throttling valves at the 
turbine and compressor inlets. 

weebest COMpPressor. 

The drive turbine is a dual-flow axial air turbine with 
90% reaction. The geometry is given in Table 1. The profile 
Shapes of the turbine rotor and of the stator blades are 
identical and the blades are of constant section along the 
radius as shown in Fig. 2. The stator has 31 blades while 
the rotor has 32 (to avoid resonant excitation from wake inter- 
ference). The two parallel stages of the turbine are designed 
for the following output and total inlet conditions: 

Pressure Ratio: 2.8 

Total Inlet Temperature: 640°R 

Flow rate: 10.85 LBM/SEC 


Horsepower: 485 HP 
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The compressor presently under test is a transonic single 
stage, axial flow compressor. It is instrumented for measure- 
ments of torque, mass flow rate, stagnation temperatures and 
pressures, case and hub wall pressures, and for unsteady 
pressure measurements in the flow field and at the walls. 

The Air Supply System incorporates an electric motor- 
driven multi-stage axial flow compressor manufactured by 
Allis-Chalmers. It can presently supply up to 12 lbs/sec of 
air at 3 atmospheres, at temperatures between 560°R and 660°R. 
The compressor is rated at 1250 HP and has a controlled 


variable speed drive. 


B. STATEMENT OF THE TASK 

The Transonic Compressor Test Rig was designed to pro- 
vide the means for obtaining experimental data in fundamental 
compressor phenomena. Following the present experiments, an 
experiment to investigate the onset of supersonic unstalled 
blade flutter is planned which would involve replacing at 
least the present compressor rotor by a rotating cascade of 
flat-plate blades. Such a rotor would not be able to produce 
the pressure ratios required to pump the required flow rates 
through the system. Therefore, it has been proposed, that 
a turbocharger compressor be fitted in series with the 
rotating cascade to provide the required flow through it. 
The turbocharger would also be driven using air from the 


Allis-Charmers air supply system. 
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In order to evaluate the feasability of the turbocharger 
installation, it is necessary to determine the mass flow rate 
required by the drive turbine to drive the test compressor 
at a given power and speed. The remaining air to drive the 
turbocharger turbine is then known and the selection of a 
commercially available turbocharger suitable for this 
application can be made. 

Thus, the performance of the air drive turbine must be 
known over the complete speed range. Of particular impor- 
tance, are the required mass flow rates for given values of 
horsepower. The problem, therefore, is to obtain the turbine 


performance map for all pressure ratios and speeds. 
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It. APPROACH 


A. BACKGROUND 

A search of the most recent literature revealed a number 
of analytical methods for the calculation of turbine off- 
Secten performance. The majority of these used in a finite 
element approach but little information on the relative suc- 
cess of these methods in practice was available. Two alternate 
methods, both used at the Turbopropulsion Laboratory, were 
those of M. H. Vavra and E. Macchi. Each was examined in 
detail. 

The method of Vavra, given in Ref. [1] is a one-dimen- 
Sional (meanline) approach using mathematical modelling and 
experimental data to express flow angles and losses. It is 
primarily a method to design turbine blading but may also be 
used to predict turbine performance for a given set of gas 
inlet and operating conditions when the blading geometries 
mee specitied. It is assumed that the axial velocity is 
constant along the blading from hub to tip. Vavra states 
that this assumtion is reasonable for blading in which the 
tip-to-hub ratio is equal to or less than 1.15. The ratio 
momileosic and 1.424 for the drive turbine stator and rotor 
blading respectively. It was thought therefore, that the 


method of Macchi might yield more accurate predictions. 
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Macchi's method is given in Ref. [2]. The method, imple- 
mented by Macchi in a computer program written for the IBM 360, 
was an extension of the work done by R. Eckert [Ref. 3] and 
Remarrison (Ref. 4]. Eckert wrote a program, following a 
Simplified three-dimensional analysis, which could be used to 
predict the performance of a single-stage axial flow turbine. 
Harrison improved the program by modifying the analysis to 
take into account streamline curvature. Both programs were 
based on the three-dimensional method developed by Vavra in 
Ref. [5]. Macchi's principle improvements to the program were 
to introduce the choice of various methods to calculate gas 
outlet angles and loss coefficients. Two methods of calcula- 
bing gas outlet angles are included; those of Ainley and 
Mathieson [Ref. 6] and Traupel [Ref. 7]. Five methods for 
calculating the loss coefficients can be selected; those due 
to Ainley and Mathieson [Ref. 6], Dunham and Came [Ref. 8], 
Balje [Ref. 9], Lonherr and Carter [Ref. 10] and Traupel 
fiet. 7). 

Macchi's computer program, as documented in Ref. ]2], was 
selected for performance predictions of the drive turbine. 

It should be noted that no card deck of the program was avail- 
able, and no results of using the program were available 


other than those included in Ref. [2]. 


Beanery ols 
The method requires the following assumptions; 
1. There are an infinite number of blades in each blade 


row so that blades downstream do not affect upstream conditions. 
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Z. the flow is axisymmetric at locations where the 
equation of motion is solved. 

Eeeeeene tlow 1s steady and ad@abatic. Thus, the total 
enthalpy through the stator remains constant along a stream- 
line and the relative total enthalpy through the rotor remains 
constant along a streamline. 

4. All equations are solved at between blade row loca- 
tions. Increases in entropy occur in the blade row upstream 
of the stations where equations are solved and the entropy 
change along a streamline between blade rows is Zero. 

>. Ihe boundary layers on the turbine casing are not 
aecounted for. 

The method of solution is as follows: 

1. Assume initial radial positions of the streamlines. 

Jeo eeaan tne axial velocaty distribution by solving 
mie caowation of motion at the stator outlet. The velocity 
distribution into the stator is assumed to be axial, and 
uniform 

BS. “Obtain stator loss coefficients. 

mM Check overall continuity and adjust the inlet Mach 
number as necessary. 

5S. Check the between-streamline continuity, and adjust 
Sereambine radial positions as necessary. 

6. Repeat this process for the rotor. 

7. Re-cycle all the above calculations, accounting for 


streamline curvature, and repeat until convergence is reached. 
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oe Et HOD OF SOLUTION 

The computer code written by Macchi was originally run 
on the IBM 360 computer. The program consisted of a deck of 
over 2000 program cards plus over 60 data cards. Since the 
deck could not be located, it was necessary to re-type the 
program from the listing in Macchi's paper. However, since 
the IBM 560 computer was soon to be replaced in the period 
in which the work was to be carried out, an alternate 
computer was sought. 

The HP-1000 series mini-computer located at TPL was 
selected for two reasons. First, the machine used FORTRAN 
as did Macchi's program. Secondly, it would be a benefit 
to TPL to have the program immediately available on the 
laboratory computer. 

The first steps were to analyze Macchi's program, in 
detail, and then to run it using his example input/output. 
In analyzing the program it became obvious that the computer 
program listing given in Ref. 2, was not the one used to 
obtain the listed output. Numerous discrepancies were found 
in the listing, some of which would have prevented the pro- 
gram from running; otners would have caused incorrect results 
to be obtained. A listing of these discrepancies is con- 
tained in Appendix E. When the program was understood and 
flow-charted, it was keyed-in at the HP-1000 computer ter- 


minal. However, modifications were required to accomodate 


Ihe) 








the program within the mini-computer disc-based operating 


system. 


D. MODIFICATION TO THE COMPUTER CODE 

Since there was no card reader, variable input data such 
as turbine speed had to be entered using data or specification 
Seaecments. This contributed in part to the most difficult 
problem, that of program size. The HP-1000 mini-computer uses 
a disc with a storage capability of 19.5 mega-bytes. However, 
the machine memory is only 124 K Bytes, of which only 29 K 
Bytes 1S available to a programmer. Also, the available 
memory is divided up, or partitioned into two 18 K and one 
ll K partitions, so that no single program can exceed 18 K. 
It was estimated that Macchi's program was over 100 K. So 
1t was clear that the program would have to be modified if 
it were to run on the mini-computer. 

The first modification was to remove all subroutines 
from the program that were not actually used. It will be 
recalled that Macchi's program contained five methods for 
calculating loss coefficients and two methods for calculating 
gas outlet angles. It was decided that only the Traupel 
method of calculating loss coefficients would be retained. 
Traupel was selected for two reasons. Firstly, it was the 
method used by Macchi in his example calculations and there- 
fore the modifed program should still reproduce Macchi's 


results. Secondly, the method of Traupel is widely respected. 
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The method of calculating gas outlet angles was totally 
Changed. Neither Ainley and Mathieson [Ref. 6] nor Traupel 
[Ref. 7] was used. Both methods required prohibitively 
large sections of computer code. The method selected was 
mmat of Vavra [Ref. 1]. 

Use of Vavra's method greatly simplified the program 
because this method predicts gas outlet angles independently 
of the inlet Mach number. Macchi's approach was to use 
Traupel's method which is dependent on the Mach number of 
the flow into the blade. 

The above simplifications reduced the program size from 
2257 lines to less than 1800 lines. However, this was still 
too large and the program could not be loaded without over- 
flowing the memory. 

The solution to the problem was found in program 
segmentation. In this process, the computer code is divided 
into a main program and several segments. Each segment is 
a "piece" of the original program. The segments are indi- 
vidually compiled and loaded. However, the segments are 
placed into memory only as they are needed to execute the 
overall program. Thus, a very large program can be made 
to run in the available 18 K partition. Since the present 
program was not originally intended for a mini-computer, 
segmentation was not straight forward. The method finally 


arrived at is detailed in Appendix C. Basically, the main 


a 





program consists of all the subroutines, while the three 
other segments contain coding which enables program flow to 
proceed in a logical manner. 

Successfully segmented, the program was run using Macchi's 
input. An output was obtained which agreed almost exactly 
With Macchi's results. All output quantities were within 1% 
of Macchi's quantities. The differences were, in all proba- 
pitty, due to the different method of calculating gas outlet 
angles. 

After verifying Macchi's program, the drive turbine geo- 
metry was input and the program was run for a given set of 
operating conditions. The results are discussed in the follow- 
ing section. Note: The "verification" of Macchi's program 
amounted to verifying that the computer code now loaded into 
the HP-1000, was indeed Macchi's code. It was not known 
whether Macchi's output data were a good or bad prediction 


of performance since they were not compared with test results. 
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ITI. RESULTS OF AXIAL TURBINE PREDICTIONS 


A. USING BOTH COMPLETE AND MODIFIED PROGRAMS 

The drive turbine geometry was input and the following 
solution flow path was selected: 

1. Stator and rotor loss coefficients were functions of 
pressure ratio. 

Eine plockage factor, &*, used in the equation of con- 
tinuity was equal to the total loss coefficient. 

Four operating points were selected to test the validity 
of the program. Three were off-design points at which measured 
data were available and the fourth was. the design point itself. 
Table II contains details of the selected test points for 
Run 1. 

The program variables were then changed and the following 
new solution flow path was selected: 

me otator ana rotor loss coefficients were those Catcu- 
lated by Traupel's method. 

Peeamie plockage factor, &*, was equal to the profile loss 
coefficient. 

After reviewing the results of Runs 1 and 2, a further 
modification was made to the program. The original program 
Contained a subroutine which checked between-streamline con- 
tinuity. If the total mass flowrate at the stator and rotor 


exits was not evenly divided between the five streamlines, 
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the radial positions of the streamlines were adjusted and all 
steps were recalculated using the new streamline positions. 
Hence, for Run 3, a subroutine was removed and the main pro- 
gram was modified so that between-streamline continuity was 


not examined. 


B. COMPARISON WITH MEASURED DATA 

ieecemesults of Run 1, 2, and 3 are tabulated in Table III. 
Run 1 showed predictions of mass flowrate which departed about 
6% from the measured data. However, the horsepower predictions 
were off by as much as 16.17%. Furthermore, the computer pro- 
gram was unable to reach a solution for the design point. 

Run 2 produced worse results as is evident from the table. 
Again, the program was unable to converge to a solution at 
fmgemadesign point. 

Run 3 produced more acceptable data. Additionally, con- 
vergence to a solution was noticeably faster and a solution 
was obtained at the design point. Because of this, the method 
used in Run 3 was used to map the drive turbine performance. 
The computer program used to obtain the results of Run 3 is 
described in detail in Appendix A and is listed in Appendix G. 
The results of Run 3 are shown plotted in Figures 3 through 8. 

To obtain the plots in Figures 3 and 6, a value of the 
total inlet temperature was approximated by the method of 
Vavra as contained in Ref. [14]. It was assumed that the 


Static turbine discharge temperature should not be less than 
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45°F (505°R). This corresponds to the approximate temperature 
at which condensation of moisture in the air, assuming 100% re- 


lative humidity, will occur. The inlet temperature was given by 


Total Inlet Temperature = DEercmOuUC lous Nempemaulre 
a a; 

where n., the total-static turbine efficiency was assumed 
eorpe S1%, and Sp» Was the total to static pressure ratio. The 
total inlet temperature corresponding to each pressure ratio 
1s given in Table I[V. 

The computer output corresponding to each point on Figures 
3 through 8 is contained in Appendix F. Only one side of the 
dual flow turbine was analyzed, thus, the resulting printed 
values of horsepower, referred horsepower, moment, referred 
moment, flow rate and referred flow rate must be doubled to 
obtain the actual turbine characteristics which have been 


plotted in Figures 3 through 8. 
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LV .. SB ESGUSSION 


mace agreement of both the predicted flow rateand the horse- 
power obtained in Run 3 with turbine test data was encouraging. 
It is to be noted however, that this agreement was obtained 
using a procedure which was conceptually incorrect. In Runs 
1 and 2, between-streamline continuity was checked and the 
Streamlines were adjusted as necessary. In Run 3, between- 
Streamline continuity was not checked, and as a result, the 
mass flow rate between streamlines was not precisely 25% of 
the total flowrate. It is noted however, that the deviations 
were less than 10.0% and while the radial positions of the 
Streamlines varied by 10.%, the differences between predicted 
and measured output horsepower decreased from 24% to 4.5%. 
Since the enthalpy change on each streamline was computed 
using Euler's turbine equation, the total horsepower obtained 
Eyeimtegration is sensitive to the streamline radial positions. 
On the other hand, the calculation of the overall mass flow- 
rate is primarily a function of the blade throat openings and 
inlet conditions of the flow. Consequently, in relaxing the 
requirement for between-streamline- continuity, the output 
horsepower was changed significantly, while the overall 


flow rate was not. 
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Using this procedure, which preserves overall continuity, 
a performance map for the turbine was produced (Fig. 3-8) 
which agreed well with the off-design performance measurements 
made at lower speeds (Table III). It is noted however, that 
the inability of the program in its original form to predict 
the measured turbine performance is not explained, and both 
the program itself and the data input for the geometry should 
be closely re-examined. 

The difficulty in obtaining convergence to a solution at 
some Operating points above the pressure ratio of 2.0 is 
meety tO be the result of choking occurring on one Or more 
of the streamlines. This was suspected but not fully explored. 

Finally, although the program was eventually made to run on 
the mini-computer, the time required to put the program into its 
final form was excessive since the original program was not 
written with segmentation in mind. When the segmented pro- 
gram was completed, only one operating point per run could be 
obtained. Thus, excessive time was spent compiling and loading 
the program. The execution time for the program averaged 2 
minutes at the lower vressure ratios and up to 30 minutes at 
the higher ones. This would be unacceptable if many points 


were to be examined. 
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V. CONCLUSIONS 


The program for calculating the performance of a single 
stage axial turbine reported by Macchi was revised, corrected 
and segmented and made to run on the Laboratory mini-computer. 
When applied to the geometry of the air-drive turbine of the 
Semenessor test rig, selecting specific options for the repre- 
sentation of loss coefficients, the revised program failed to 
converge when design-point test conditions where input. Also, 
the computed horsepower was in error by as much as 24% when 
the program predictions were compared with specific test data 
obtained from the rig at off-design (lower speed) conditions. 
The revised program did however closely reproduce the results 
paven by Macchi in his original report for a specific turbine 
geometry. 

When the requirement that the computed stream surfaces 
be such that they divided the flow exactly into equal 25% 
increments was removed, the program converged satisfactorily 
for design point conditions and gave agreement with test data 
to within 5% in flow rate and horesepower at off-design 
Soma tions. 

The complete performance map for the air drive turbine 
was obtained with the program following this revision. Based 
on the favorable comparison with data so far obtained, the 
map is likely to describe the performance to better than a 10% 


uncertainty. This is considered to be satisfactory for 
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sizing the turbocharger for the proposed compressor rig 
modification. 

The following recommendations are made concerning further 

application or development of the computer program: 

1. The failure of the program to converge before the 
final revision was made should be analysed closely, 
and the final revision removed if possible. 

2. The geometrical input for the air drive turbine (which 
was taken from drawings) should be reexamined and the 
physical dimensions of the blade rows themselves should 
be measured. 

ee Gensideratiom should be given to putting the corrected 
original version of the program onto the IBM 370 com- 
puter so that, when successfully operating, a turbine 


map can be calculated with a single load. 
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TABLE I 


TURBINE GEOMETRY 


(see Figure 2; Dimensions in inches) 


























STATOR: 

Hub Radius ss—<SstsSsS——Ci TH 
Mean Radius 5. 196° 

Tip Radius 5. Our 

Blade Chord | 100s 

Blade Suction Side Radius of Curvature 2.8065 
Maximum Blade Thickness ot Oe 

T.£. Projected Thickness 20'S 

T.E. Normal Thickness .0186 

ROTOR: 

Tip Radius Siac on 

Blade Chord LCOS 

Blade ‘Suction Side Radius of Curvature Zo 005 
Maximum Blade Thickness SD 

ieee cojyected Thickness 205 
iteeeenermal Thickness .0186 

Tip Clearance .Ol(estimated) 
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TABLE III 


COMPARISON OF PREDICTED TURBINE PERFORMANCE 
VS MEASURED PERFORMANCE 

















RUN I 
FLOWRATE HORSEPOWER 
POUNT Jer leche 10141 Ot ae MEAS. *+DIFF. FRED eG lye MEAS DIFF. 
1 00 9.942 6.09 OOS aA Oise or 
Z 4.74 4.698 Oms9 2.0 OS Oe / 
5 7.04 (OSS 0.009 163.64 ace) 4.86 
4 Ne Gs 10.85 ae Nees 485 Yeats 
RUN 2 
i 6.06 Dag oS a S032 BLO eel | 
Z 4.90 4.698 4.29 49.76 O522 7) Pca 
5 50 VO Saou sO %7 0 1720 Zo oT 
- NGa LO: 8.5 ee Nace 485 Sass 
RUN 3 
1 oe Date Se Lt Srlyz Ge 2.74 
2 4.66 4.698 Ook Gio OS 220, 2.09 
3 7.04 (eS S 0-10 E/ 95.08 ae 4.47 
+ 10.40 Oe SS 4.15 444.18 485 8.42 





vas 


C: Computer program would not converge to a solution after 
a large number of iterations. 








DABLE. Ly 


VALUES OF ASSUMED TOTAL INLET TEMPERATURE FOR EACH 
PRESSURE RATIO. GIVEN IN FIGS. 3, 5, 6, AND 7 


PRESSURE RATIO LOLA INCE TEME ERA Ube CR) 
Lod 945.55 
eo 027-0 
| es! Bl Tac 
Z 0 591.0 
Zc 0052.6 
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a5 6 b216;:,0 
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CY Boz E-L ATM. 3 ATM. 
To: F30°R }O0-200 °F 










AIR DRIVE 
TURBINE 


Z 


EXISTING 
SYSTEM 








~_—E— == — eS ae eee Eee a eee 


PROPOSED 
MODIFICATION 


SNe KR SO- 
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Og 


Peak i: SCHEMATIC.<CF THE COMPRESSOR TEST RIG, WITH 


PROPOSED MODIFICATIONS 
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APPENDIX: A 


DESCRIPTION OF THE COMPUTER PROGRAM 


A-1. INTRODUCTION 

To enable the program to run on the laboratory computer, 
the program was divided into 4 parts; a main program and 3 
segments. A detailed discussion of program segmentation on 
the HP-1000 computer series is contained in Appendix C. In 
the description which follows, the program is treated as if 
it were one large program with many subroutines. 

The description follows the individual steps from 
start to finish in the analysis. A program flowchart is 
given in Figure A-1 and the FORTRAN symbols used in the 
program are listed in Tables A-I to A-IX. 
h-2. DESCRIPTION 

rae) input Data 

There are 4 basic categories of input data; tur- 

bine geometry, operating conditions, special data and pro- 
gram control parameters. Since there was no card reader 
input device on the computer, all datawere entered using 
either data or specification statements. Explanations of 
the turbine geometry, operating conditions, special data 
and program control parameters are found in Table A-I through 
A-V. The nomenclature for the blading is given in Figure 
A-2. 
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mead initial Geometric Calculations 

The first calculation performed is to establish 
mens streamline locations at the stator inlet (station 0). 
The streamlines are initially positioned such that there are 
equal areas (25% of the total flow area) between them. Next, 
blade heights of the stator and rotor are calculated using 
the hub and tip radii of each blade. Blade spacings for the 
stator and rotor are computed at 3 streamlines; hub, mean 
and tip. The blade spacing on the mean streamline for the 


Stator is given by 


eee Ti 
S = ph (A-1) 
where S = Blade spacing 
Z. = Number of stator blades 


Rm = Mean stator radius 
A-2.3 Calculation of Gas Outlet Angles 
Subroutine VAVRA calculates gas outlet angles 
BOEeoetn stator and rotor. The method is that of M.H. Vavra 


[Ref. 1]. The equation programmed in the subroutine is 
-i f @ 
=I OW +4 te (| — £28 a 
oO = Coit S S GO (A-2) 


Gas outlet angle 





= 
ae 
M 
ry 
@ 
{2 
iT) 


m 
il 


Throat opening 
S = Blade spacing 
t,= Projected trailing edge thickness 


This method is much simpler than that used by Macchi since 
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there is no variation in outlet angle with Mach number (for 
sub-sonic conditions). Therefore, once calculated, the stator 
and rotor exit angles remain unchanged. Subroutine VAVRA 
computes exit angles for the hub, mean and tip streamlines. 
The outlet angles at streamlines two and four are computed 
later in the subroutines STATR and ROTO2. 

Before printing the input data, the program 
calculates the mean throat opening for the stator and for the 
rotor. The ten equally spaced radii and corresponding throat 
openings (part of the input geometry) are fitted with a fourth 
order Ghebyschev polynomial. A throat opening corresponding 
to the mean radius is thus obtained. In the present appli- 
cation of the program to the drive turbine, the mean throat 
opening was obtained from the design drawing of the blading 
Shown in Figure 2. It was assumed that the throat opening 
varied linearly with radial position and hence the throat 
openings at other radii could be calculated. The resulting 
throat openings are shown in the computer output under the 
heading of "Input Prints". The design values of the stator 
and rotor throat areas were obtained from the original 
design notes of M.H. Vavra. 

A-2.4 Calculation of the FlowRate 

Subroutine CHAN is called to calculate the mass 

flow rate entering the stator. The equations used are as 


follows: 
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Tro 
| + ee Per 


V = 4 GceORT a) 





(A- 3) 





P 
be sea (A-5) 
ee Me 
p= P/RT me 
2 2. 
A= 7 | Rip am 108 | ae 
ae (A-8) 
Meee = oo | Rie (A-9) 
: TO Cc 
Mot is the reference (dimensionless) flowrate 


fiers used to check overall continuity later in the progran. 


A-Z.5 Solution of the Equation of Motion for the Stator 


Subroutine STATR is called to solve the equation 
of motion for the stator outlet conditions. The equation of 


motion which is programmed is as follows: 


AE) poeta] (warm $2) (HeeSa%) 





AX i ie 
tor | ee ee 2 ty SS 
YG 3) TAN XY ——— dX | SIN OX 
Ci Cos7el aH Ci H cos‘ol 2 ENE 
| or - SS So a 
fae Vicu a aA‘ t Yi Var, AX,  (A-10) 
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where Gi = Ae u (a constant to convert H, the enthalpy 
from BTU TO bili 
LBM sec 


“¢ c Vane) 


VV | e Amat velocity wit a streamline 
a. us Axial velocity at mean streamline 


x - ROT) . otreamline radius 
Fag Mean streamline radius 
Y, = Vee 
ds” _ a |p | 2 cit costa 
es Ol Pe 2 > 
aX, al X | - Yi * Vau 


Ci H cos*x (1~ @ ) 
Ge soit aeein Loss (COC NETene 
(which is initially assigned an 
estimated value) 
The derivation of this form of the equation of motion is 
given in Appendix B. However, at this stage of the analysis, 


the streamline curvature is assumed to be zero. Therefore, 


the equation of motion becomes: 


_d (dn %") 29 Tavielp 2c . 


a eC Cos, 
— SIN & + —S— 
d Xt AX Kt ? Y7 Vai 


dH A oe C,H COs" ds* | 
aX, | Giese Ver m AX oe 


The equation of motion is solved when the value 

of Y at each streamline satisfies the equation. The solu- 

1 
tion is to first put the equation in the form: 
2 
Yi 
d (Ln ) = TL (x) 
AX, 


where I(X) consists of the right hand side of equation (A-11). 


(A-12) 


Integrating equation (A-12) yields; 
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X 


may * = { WSO ai +A C (xe 13) 


XO 
where &£nc* is the constant of integration when x = 1 and Y-= 
1 


ietth these boundary conditions Eq. (A-13) gives 


Ln c* = i> ( Lixd xX: (A-14) 
XO 
Meme Eq. (A-14) in Eq. (A-13), 


ay *. 


Xo 


( | 
IO)dX-| LOdd X (A-15) 


which becomes m 
{ 


a s 
Un, = ( L(x] dX oe 
( 
Taking the inverse natural log and the square root of both 


; x 
side . (*rudx 
af, 
oS (hai 


Equation (A-17) is the form of the equation of motion solved 
in subroutine STATR. Solution of the equation gives five 
Values of yy Fis eMcmencavalveaG «the idealal velocity at 
each of the five streamlines. Initially, the value of Y, 

1s taken to be 1 and the value of a is taken to be zero. 
In succeeding iterations, the eames value of Y 1s used 
to obtain a new value of —: and so on. 


1 
After calculating five values of Y> Chem sea cor 


exlt conditions are calculated at each streamline from the 


geometry of the velocity diagram. The convention for positive 
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and negative angles and velocities is defined in Figure A-3. 


The required relations are the following: 


Va, = Va, (3)- Yi on 
Vin, = Va, > TANS, (A-19) 
V, = Va J 60s & (A-20) 
Ve, = - Va, | AR /av | (A-21) 


where L is the axial distance between stations and AR is the 
change in radial position of the streamline. Vee the radial 


component of velocity, is taken to be zero at this stage in 


| vaee Vie 


the calculation. 


V = 
V a (A-22) 
a = T +o - + 
AGe JCp (A-23) 
ans = ie. _ Tro-T1 
ie z (A-24) 
PR. = P, ; Pro a CxS 
p Es P Ti r-{ 
a a Tro (A- 26) 


M, = V / | Vac RT (A-27) 


After the above quantities have been calculated at each 


streamline, subroutine STATR returns to the main program. 
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A-2.6 Calculation of the Stator Loss Coefficients 
ihe calculation of the staton) loss coefficients 
at each streamline is accomplished by subroutine ALOS1. 
The method of solution to obtain these loss co- 
efficients is that formulated by Traupel [Ref. 7]. I50) 
ieaupel s method, the value of the total loss coefficient is 


given by 


=e cae teas (23) 
total profile wall remaining 


mime calculation of & requires 9 subroutines. Figure A-4 


total 
describes the connection between the subroutines and subroutine 
EOS! . 

The first step is to obtain the value of the 
eecal profile loss coefficient, a ED 1s defined by Traupel 


to be 


5, ‘ Fro Xm X § i to a ¥ (A-29) 


where EDO = initial value of the profile loss coefficient 


ie mach number correction factor 
a trailing edge thickness correction factor 
Em = loss coefficient due to mixing losses and separa- 


tion losses 
Ef = loss coefficient due to fan losses 
ties totalmprotile loss coefficient is calculated 
in the following manner. First, data for initial profile loss 
(E) as a function of gas outlet angle (a,;) for various values 


of gas inlet angle (a,) is read from an array (Fig. A-5). 
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is 15 done by subroutine TRAU1] and functions XPO and YC. 
The values of EDO are contained in two arrays XPO1 (5, 8) and 
XPO2 (6, 8). This is because the data shown plotted in Fig. 
A-5 has been divided into two sets. One set is for values 

of a, between 40° and 80°. The other is for values of a, 
between 80° and 170°. The FORTRAN symbols for the two ranges 
of values of a, are ALFO1(I) and ALFO2(I) respectively. The 
FORTRAN symbol for the gas inlet angle is ALFl (J) once the 
data pints selected from the plot are entered, fifth and sixth 
degree polynomials respectively are fitted through the data 
points. The value of EDO can then be determined for given 
values of a; and ay. 

Pre nach number correction, a is obtained from 
Bae. A-5. Subroutine CSIM calculates the value of = using 
eeraighnt line approximations of the plot. 

Subroutine CID calculates the remaining terms in 
the expression for ED’ These are Xo em? Ee: They are ob- 
tained from the data in Fig. A-6 using the linear interpolation. 
The abscissa of the curves for i and em Se) CanGh Che fe  Oieet eae 
where f is defined as S 


ET SINK, (A-30) 


where 6 Noumaletraiding edge thickness. 


blade spacing. 


ct 
i 


a; = gas outlet angle. 


iMeomiloss. coetticient due to wall friction, ey? 


is calculated using 
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¢ = fo Xe Tb SINA 


L (A-31) 





where t - blade maximum thickness 
£ = blade height 
This equation is programmed in subroutine CSIW. 
The value of ER 1s obtained using subroutine CSIR. 
ER is defined by Traupel to be an all-inclusive loss coeffi- 
cient which accounts for any remaining losses not previously 


Metined. it is written as 


f = XX. f (A-32) 


ERO MS an initial value of ER which depends on the value of 


@?, where $¢ is given by 


Vi SINS, 
D 7 V (A- 33) 


in which v, = true velocity of gas 


V blade speed 
Ee plot of ERO vs 6 1S shown in Fig. A-7. The correction XxX, 
1s a function of s/2 where 

s = chord length 


Q 


blade height 

PigelssObtained using the data in the lower half of Fig. A-7. 
The total stator loss coefficient 1s computed 

for 3 streamlines; those at the hub, mean and tip. 

The loss coefficients at streamlines 2 and 4 are 


obtained by linear interpolation. 


>i 











A refinement to the stator loss coefficient may 
be applied depending on the input value of one program control 


bamameter. The following 3 variations of b . are available: 











Piel 
+ §, 
¢ _ Peer _ | 
a — | os 
r 
2 
Fro (A-34) 
¥ : &. (A-35) 
and | ol 
r 
| + f., | 
|+ §o8* 
> 7 ) ee 
sey AC (A-36) 


where € = loss coefficient calculated using the method of 
0 


Traupel y-! 
o~ 
: (9) 
eye lt OR ae 


(A-37) 


The values of the program control paramenter 


required to select between options are given in Table A-V. 





Before returning to the main program, subroutine 
ALOS1 calculates a value of €* which is a blockage factor to 
be used in the equation of continuity. There are three ways 


to define €*; they are as follows: 


¢ ; = t (A- 38) 


i. 4 
“ f. (A=39) 
$ 


f i. t, (A-40) 


Peon ee ooOluicion of the Continuity Equation After Re- 
turning to the Main Program 
ime overall Continuity at the stator exit is 
checked. Subroutine FLOWR performs this task. The flow chart 
POreGOWR 2S given in Fig, A-8. In FLOWR the mass flow rate 
required by continuity is checked against the calculated mass 
flow rate. If the calculated flow rate does not agree with 
that required by continuity, adjustments are made to the axial 
velocity and/or the inlet Mach number, as will be explained. 


The mass flow required by continuity is 


Wire M REF 


Mesan 2 
Zs Am K m (A-41) 


where Meer = reference mass flowrate as computed in subroutine 


CHAN 


ao 








Meath oOrt Stator blades 


S 
i = mean stator throat opening 
R, = mean stator radius 


The mass flow rate at each streamline com- 


puted in this subroutine is 





a Pre Tte ACT) 
Maer | P T, z 
where Z is an area reduction coefficient defined by 
KK 
>. S 
¥RK %. 
H Sealed f eas 


Z gives the percentage of flow area hetween the blades over- 
which it is permissable to assume a uniform velocity. The 
boundary layer on both sides of the flow limits the available 
flow area and the backage facotr, Z. accounts for this. 
Equation A-43, Z is seen to be afunction of the energy para- 
meter H*** and &*. &* is the value of the loss coefficient 
returned from subroutine ALOS1. The energy parameter is 


defined as 
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Oo 
Cn 


Energy thickness 
KKK = = = 
sa Displacement thickness CA 8) 


H*** can be written as 


O 
}— 




















| 4 Y 2 Y 3 Y “ 
p P Een eR eee 
Pee Xer-l Smt Sm-+\ Tm +| FT m+! Um +) 
2: 3 y 
} | 
Yer! m+! 3m+tl Smt 7 mri Tmt 
(A-45) 
where: 
m= .15 
x = ere — for unchoked flow 
E TO 
Vek. 
x mele aeeet iu for choked flow 
and ; Y_ 
ii ieee 
Porit ~ | yi! 


The derivation of Z and H*** os given in Appendix B. 
The expression for %, the flow function, 


for unchoked flow is 


(A-46) 





and for choked flow is 
OD = 2 2t 
rel y+ 


SS 





(A-47) 





After calculating for each streamline, the flow 
rate 1s integrated from hub to tip and the resulting value is 


compared with m If the two values of flow rate agree to 


reqd° 
within a specified tolerance (see Table A-IV) continuity is 
Sememaered to be satisfied. Then, after calculating the total 
percentage of mass flow between adjacent streamlines, sub- 
routine FLOWR returns to the main program. 

If the flow rates are not within tolerance the 
program checks to see if the actual mass flow is to high. If 
it is to high, the value of the axial velocity is lowered 
proportionally to the difference between the actual and 
required flow rates. 

If the actual flow rate is too low, the procedure 
is more complicated. First, the flow is checked to determine 
whether choking has occurred. Streamlines one and five are 
checked. If the flow is in fact choked at those streamlines, 
the inlet Mach number is lowered and the program loops back 
to recompute the reference mass flow rate and repeat the 
complete procedure. 

If the flow is not choked, the axial velocity is 
raised proportionally to the difference between actual and 
required flow rates and subroutine FLOWR returns to the main 
program. 

A-2.8 Calculation of the Rotor Inlet Conditions 

Continuity having been satisfied through the 
stator, the rotor relative inlet conditions are calculated. 


In subroutine ROTO], the following expressions are used: 
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U, = WwW ; hiro. 
[a R RoTOR 
War, = Va, - U 
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6, = Wew [Men] 
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W, = Ve. : Ww 
od 
7 e = ( Ty ee AY, } 
' ee | 
Aqe dep 
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ee — ( Tre )(-24) 


Where Trp» POE and He are equivalent temperature, pressure 


(A-56) 


and enthalpy respectively. 
noe Calculation of the Rotor Exit Conditions 
Calculation of the rotor exit properties follows 
the same procedure as was used to compute the stator exit 
properties. The process is outlined here with notable 
Gitterences explained. Subroutine ROTO2 calculates the 
rotor exit properties. A flowchart of ROTO2 is given in 
Eee A-9. 
Mic ttarstestep in ROPOZ is to solve the equation 
of motion for each streamline. The equation of motion in 


terms of relative quantities is 


A(dnYy") ser = GE) 
dX 2 gre ; 


ds* aba Ea Z 4 Um CoS 
ay | A TAN a - 2 sna - = ee Ca. 





BD ORew> cos” bez G, cos Ba ; od He 
Wace Va, Yo Vay qd X2 


| - Ci He cos’ Ba . ds." (A-57) 
Vai Vas ol Xz. 
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At this point in the calculation, streamline curvature 1s 
Meetected. Hence, Eq. (A-57) reduces to 


a (Ln Ya) ol 62 a 


= eae eee 2 a Se a == 
AX dTran Bo 1X, x Sitwoo 


4 Um cos Ba Sin By _ 2 Um Uz cos Ba 
Ya Vay Ya Vary ° 


Cicos"Br ofHe e- C. He cos762| Sz 
x." Ves dX2 Wee Vax ONG: 


(A-58) 


The derivation of Eq. (A-57) is contained in Appendix B. 
Equation (A-55) is similar in form to Eq. (A-10). Hence, the 
method of solution is identical to that employed by subroutine 
STATR. However, after solving the equation, the value of Y, 
at each streamline is examined to determine whether or not it 
milis into the range waa 60: Values of Y, greater than 
2.0 are set equal to 2.0 while those less Ghai. cease Ser 
pale to .2. Successive values of ty at each streamline are 
compared, and when the values of successive iterations are 
within a specified tolerance (see Table AerliN tne el tengac Lon 
ends. The values of Y, are used to calculate the rotor exit 
conditions using the following equations: 


va = Vag (3) Ts 


(50) 
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COS Ba (A-60) 
We (- Va, )'D+ CL 
5 2 (A-61) 
2 
Ta = Tre ~ bid 
Ade I &p (A-62) 


< 
é 
jo 
| 
os 
3 
Cc 
“ 
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eS) 
a= iy el) 
= we (A-64) 
Tas = Tte _ Tre - Ta 
e ee (A-65) 
ae 
re 
p Tas 
* Pre Tre (A-66) 


Subroutine ROTOZ, then returns to the main pro- 
gram. 
Mteer cameulating the rotor outlet conditions, 


the rotor loss coefficients are computed. Subroutine ALOS2 
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calculates the rotor loss coefficients following the process 
used in subroutine ALOS1 for the stator losses. The principle 
exception is that a tip clearance loss is also calculated and 
added to the total loss coefficient. The tip clearance loss 
coefficient is obtained from subroutine ALEAK which uses a 
straight line approximation to the curve shown in Fig. A-10. 
Subroutine ALOS2 also computes values of &* and one of the 
three refinements to ER° 

subroutine FLOWR is called to check continuity 
at the rotor exit. [If continuity is satisfied, the program 
continues. If not, the same procedure is followed as pre- 
viously explained for the stator outlet (Fig. A-1l). 

A-2.10 Accounting for Streamline Curvature 

All calculations to this point have neglected 
streamline curvature and assumed that the streamlines remain 
miced through the stator and rotor (Fig. A-11). The radial 
shift in a streamline between stator inlet and rotor outlet 


can be written as 


AR= eae — RK RoToR 


INLET OUTLET 


(A-67) 


This is the net radial shift in a streamline between stations 
Meumiicteateor inlet) and '2' (rotor outlet). It is shown in 
Section 16.4 of Ref. [5] that the radial shift in a streamline 
between the stator and the rotor (station 1) can be written 


as 
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[ = 
§ — K STATOR ~ 2. K stator K KotoR 
OUTLET INLET OUTLE | 


(=o) 


The angle between the meridional velocity Ne and the axial 


velocity Vy is A. The radial velocity ee can be expressed as 


Ve = Va Tan d 


ena trom Fig. 16(1) of Ref.[5], it follows that 








= AGR 
Tan d = Se 
Meane Eq. (A-68) in Eq. (A-67), 
AR 
--Ve 
VR aL 
where AR . Average streamline slope 


ZL 
melso, from using Eq. (A-68) 


ot L 





os \ = 
a Jar? (QL)” 


-Rearranging; 





(A- 69) 


(A-70) 


(A-71) 


One i a) 


Ne7 3) 





The remaining term used in the calculation of streamline cur- 


vature (Section 16-4 of Ref. [5]) is 


6R 
aa 


where K is the so called curvature factor. It usually has 
a value between 4 and 6 and in the program its value is taken 
to be 5. Having calculated cosi, AR and 6R, the program 
repeats the solution process. However, the only quantity 
which is unchanged is the reference mass flow rate Meee In 
subroutine STATR the equation of motion is solved, this time 
accounting for streamline curvature. The same is true in 
subroutine ROTO2. 

The flow path of the program is identical to the 
section which did not account for streamline curvature. Next, 
the program computes an average pressure ratio at the rotor 


outlet using the expression 

















Po, = Pa. 7 ae Py 
Pro Fo | STREAMLINE al Pro oe 
2 
R. Pa. Po 
es + 
Pro Pro ee 
Sb as S.L- 
3 5 s 
(A-74) 


If this pressure ratio is within a specified tolerance to 


the actual pressure ratio (which is input data) the program 


oo 





Bioecccas, tO the final stage of the calculations. [f not, the 
inlet mach number is adjusted by an amount which depends on 
@iemdicrenenee between the calculated and specified meee 
Baetvos. it the calculated pressure ratio is too high, the 


Mach number is lowered using 


z Pee sotimemna tl Oslin Terence - 
sO Mo 18 (A-75) 


If the computed pressure ratio is too low, the Mach number 


meraised wsing 


owe + Pressure Ratio Difference 


ss . (A-76) 


In both cases, the program loops back to subroutine CHAN and 
proceeds to compute a new reference mass flow rate based on 
the new value of the inlet Mach number. The entire process 
is then repeated until the pressure ratios agree within the 
eoeeiined COLerance . 
ao rinal Calculations 
Stator and rotor outlet conditions not previously 


calculated are computed as follows: 


-| Vie 


xX a TAN 
Va 





Ce) 
a 


Var 
V2 a (ars) 
COS Az 
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Va = V5 . We, 





Ak aes I Goege Ue 
qeN 
= Tro a dh 
ot 
v 


Pr, 





\\ 
4 
~ 
tana 
— 
pid 
Sa) 
i 
ee 





a. 
T-l 
Pr = Ir Tro 
Ty 
v-! 
Ais Tro fa 
Pro 
ROTOR EXIT _ Wa 


RELATIVE MACH # 
¥ Rae le 
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CNS 


=763))) 


-81) 


82) 


- 84) 








Tr, = T Pr | 
Tis ce 
lie. 
- = Tro - Tra 
Tto =i linre 
/ = Ito- TT 
T-S . 
Tors rye 
Stator Blade Tei 
Eftaerenc S 
— _ ie eels, 
iciency 
Tro -Tais 
* 
Yr = lias {21s 
Tro -Teais 


Head Coefficient = 


Blade-Jet Ratio = 


Stator Exit 
Relative Mach # 
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2 9¢ T Cte -Tais) 


fo 


all 


[Head Coefficient] 


Wi 


(A- 86) 


(A- 87) 


(A- 88) 


(A-89) 


(A-90) 


(A-91) 


(A-92) 


(A-93) 


(A-94) 





The turbine horsepower is obtained by integration. 
The Ah term at each streamline is weighted by the percentage 
of mass flow at that streamline. The product is then niece 
grated from hub to tip and result, Ah, is used in the tur- 


bine horsepower equation 
Ah] - pA 


5 So (A-95) 


la 
The moment is calculated using 


(H.P, )( 550) 


OC) (A-96) 


Referred horsepower, moment, mass flowand RPM are calculated 


using 


H.P. 
‘REF SO iG (A-97) 


MA REF = eee (A-98) 


ier = Sa a (A-99) 


imiuegee = > mee (A-100) 
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where Tro 


= i 519.4 





Pro 
aha 


ine vwalwessatcine total-static efficiency, total - 
mecal efficiency, total-static pressure ratio, total-total 
pressure ratio, head coefficient, blade/jet ratio, r* and 
inlet mach number are then averaged. 

With all calculations completed, the results are 
mrinted under the heading "STATOR SOLUTION", “ROTOR SOLUTION", 
and "OVERALL TURBINE CHARACTERISTICS". 
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TABLE A-I 


TURBINE GEOMETRIC INPUT DATA (STATOR) 
(see Figure A-2; Dimensions in inches) 


FORTRAN SYMBOL DESCRIPTION 

oS Number of blades 

RS (1) Hub radius at stator outlet 

RS (3) Mean radius at stator outlet 
RS(5) Tip radius at stator outlet 

c Blade chord (mean) 

GI Blade chord (hub) 

CO Blade chord (tip) 

E Blade curvature (mean) 

Ee) Blade curvature (hub) 

EO Blade curvature (tip) 

7 Maximum blade thickness (mean) 
TL Maximum blade thickness (hub) 
TO Maximum blade thickness (tip) 
RE Projected T.E. thickness (mean) 
TEI Projected T.E. thickness (hub) 
TEO Projected T.E. thickness (tip) 
TN Normal T.E. thickness (mean) 
TNI . Normal T.E. thickness (hub) 
TNO Normal T.E. thickness (tip) 
Al(1-10) Ten values of throat diameter at 


OMe cituadeley spaced radii 
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FORTRAN SYMBOL 
AL 

ALI 

ALO 

Re (1) 

me ( 3) 

RC(5) 


DESCRIPTION 

Blade camber line length (mean) 
Blade camber line length (hub) 
Blade camber line length (tip) 
Hub radius at stator inlet 
Mean radius at stator inlet 


Tip radius at stator inlet 
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TABLE A-II 


TURBINE GEOMETRIC INPUT DATA (ROTOR) 
(see Figure A-2; Dimensions in inches) 





FORTRAN SYMBOL DESCRIPTION 

ZR Number of blades 

RR (1) Hub radius 

RR(3) Mean radius 

RR(5) Tip radius 

CR Blade chord (mean) 

CIR Blade chord (hub) 

COR Blade chord (tip) 

ER Blade curvature (mean) 

EIR Blade curvature (hub) 

EOR Blade curvature (tip) 

TR Maximum blade thickness (mean) 
TIR Maximum blade thickness (hub) 
TOR Maximum blade thickness (tip) 
TER Projected T.E. thickness (mean) 
meiR Projected T.E. thickness (hub) 
TEOR Projected T.E. thickness (tip) 
TNR Normal T.E. thickness (mean) 
TNIR Normal T.E. thickness (hub) 
TNOR Normal T.E. thickness (tip) 

ci PC Tip clearance 


71 





FORTRAN SYMBOL 


A2(1-10) 


ALR 
ALIR 
ALOR 
CV 
CK 


DESCRIPTION 


10 values of throat diameter at 
10 equally spaced radii 


Blade camber line length (mean) 
Blade camber line length (hub) 
Blade camber line length (tip) 
Axial distance between stations 


Curvature Factor 
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TABBE A-ILt 


TURBINE OPERATING CONDITIONS (INPUT DATA) 


FORTRAN SYMBOL DESCRIPTION 
AMC Assumed inlet Mach number 
Assumed stator exeit Mach number 
AMS (absolute) 
Assumed stator exit Mach number 
AMR (relative) 
PTO Doerr ret pressure (Po) 
TTO Total inlet temperature ('To) 
PR Poems ate ressube sacLo 
RPM Cpenacine speed (RPM) 
VA1(3) Assumed axial velocity in stator 
VA2 (3) AS clic Cedxeldieve LOCl ey ell nolo hr 
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TAB ERAS TV 


OPECIAL INPUT DATA 


FORTRAN SYMBOL DESCRIPTION 
Toler mce for convergence of 
HO 1 : ae: 
equation of continuity 
TOL 2 Tolerance for between-s.L. 
continuity (not used) 
Tolerance in pressure ratio 
Bony 2 convergence 
TOL 4 Tolerance in equation of motion 
convergence 
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TABLE A-V 


PROGRAM CONTROL PARAMETERS 





FORTRAN. SYMBOL POSS] Bigeey vlue EFFECT/MEANING 
Won eSe nes itis Im Siub- 
IND i routines STATR, FLOWR, 
ROTO2 
i! No printing in the above 
ECL il Rotor is shrouded 
1 Rotor not shrouded 
EOZ 1 E = Fy 
6 —& = &— (Y Pressure Ratio) 
: o  Suei8 
ICON 1 = -S8roraz 
: & = SpROFILE 
: 6 = EroTAL 





TABLE A-VI 


FORTRAN SYMBOLS IN THE MAIN PROGRAM 














FORTRAN SYMBOLS SS IOI 
BESP g* = [1 + Le (.8) 7] 
OT Stator throat opening (hub) 
OO Stator throat opening (tip) 
OIR Rotor throat opening (hub) 
OOR Rotor throat opening (tip) 
0 Stator throat opening (mean) 
OR Rotor throat opening (mean) 
ANG2I Stator gas outlet angle (hub) 
ANG20 Stator gas outlet angle (tip) 
BETAI Rotor gas outlet angle (hub) 
BETAZ Rotor gas outlet angle (tip) 
G Grav. constant, 32.174 ec 
ea (eo Eben, BE 
EXP1 ee 
EXP 2 ¥=2 7, 
ERRE Gas constant, 53.3459 se 
EMME Molecular mass, 28.970 LBM/LB MOLE 
GAM yY, Ratio of specific heats 
ETAT Total-total efficiency 

ETAI Total-static efficiency 

ETAS Stator blade efficiency 
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FORTRAN SYMBOL 


DESCRIPTION 
































BLAR Rotor blade efficiency 

RSTAR Theoretical degree of reaction 
LOS Head coefficient 

BLAJE Blade/jet ratio 

DR1 Radial shift of steamlines 

AMW1 Stator exit relative Mach Number 
AMS1 Stator exit absolute Mach Number 
AMV 2 Rotor exit absolute Mach Number 
AMR2 Rotor exit relative Mach Number 
DE LH A 

HP Horsepower 

AMOM Moment 

THETA 9 

DEiaA 6 

HP1 Referred H.P. 

AMOM1 Referred moment 

RPM1 Referred RPM 

WLBM1 Referred mass flow rate 

ELAS Average total-static efficiency 
BETA6 Average total-total pressure ratio 
ETA6 Average total-total efficiency 
AKISS Average head coefficient 

RSTARS Average theoretical degree of reaction 











FORTRAN SYMBOLS 


TABLE A-VII 


FORTRAN SYMBOLS IN SUBROUTINE CHAN 


DESCRUP ION 





Teo? total temp. at station 9@ 
Inlet Mach number 

Pag» total pressure at station 9 
Streamline radii at station 9@ 

M, required mass flow, pAV 
Static temperature 

Velocity 

Static pressure 

0, density of air 

ere reference mass flow 
of M at each streamline 


oe 





TABLE A-VIII 


FORTRAN SYMBOLS IN SUBROUTINE STATR 











FORTRAN SYMBOL DE SCRI PTION 

ALFAI1 Stator gas outlet angle 

x Ratio of streamline radius/ mean 
radius 

AMS Mach Number at station 1 

if Static temperature 

P Static pressure 

V1 Absolute velocity 

VA1 Axial velocity 

Y Ratio of axial velocity to mean 
ataleve LOCiL LY 

S Entropy 

DSDX Entropy gradient between streamlines 

VU1 Tangential velocity 

PRAT (Total-static pressure ratio)! 

miis Tits 


DALF AR 








RSF Mean stator radius 

R e R 
DELR Stator in - rotor out 
ZETAPS § 

p 
ZETAS § 

Ss 
VRI1 Radial velocity 
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TABLE A-IX 


FORTRAN SYMBOLS IN SUBROUTINE TRAU2 

















FORTRAN SYMBOL DESCRIPTION 
CSIP Xp» correction to 
R re initial profile loss coefficient 
ANGI Gas outlet angle 
ANG2 Gas inlet angle 
Rl xX? Mach No, correction 
R3 Sp» remaining loss coefficient 
R2 5 oss coefficient due to wa 
We frre tion 
RPRO Sp» Eocaleprotize loss coctiilewent 
cL Rotor tip clearance 
YCL Tip clearance loss coefficient 
ROT Total loss woeffmeient 
i Blade spacing 
DEZ Normal trailing edge thickness 
HM Blade height 


trailing edge thickness 


CSID Xs> correction factor 


Sloss COCTficientc due to fan 
— £ losses 

SehoOss —COCrrI1Cc1enG cue to nxn g 
POI ™ .and. separation. 7 oS 
UM Tip speed 
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TABLE A-X 


FORTRAN SYMBOLS IN SUBROUTINE FLOWR 








FORTRAN SYMBOL DESGRIPTION 

PRATGR Critical pressure ratio 

PHICR ®opIT: critical flow function 
HSTAR H***, energy parameter 

oll Z, area reduction coefficient 
PHI 6, flow function (unchoked flow) 
ARAT Streamline throat DIA/mean throat DIA 
AS Mean stator throat diameter 

AR Mean rotor throat diameter 

WREQ M required to satisfy continuity 
WSUM M calculated 

WPER % of M at each streamline 
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TABLE A-XI 


FORTRAN SYMBOLS IN SUBROUTINE ROTOI1 


FORTRAN SYMBOL 














OMEG. w, Wheel speed (RAD/sec) 
p wR ator mean 

R 
2 ye rotor mean 

R 

stator mean 

WU1 Wu, see figure A-3 
BETAI 8,, see figure A-3 
Wl Wi, see figure A-3 
is Equivalent temperature 
Pale E@Utiarenespnuessune 
lade Equivalent enthalpy 
ZETAR 3p rotor loss coefficient 
ZETAPR fa, TOtor pLOLtale Less Cocti1elent 
DHEDX Enthalpy gradient between streamlines 
DSDX Entropy gradient between streamlines 








DESCRIPTION 
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TABLE A-XII 


FORTRAN SYMBOLS IN SUBROUTINE ROTO2 





FORTRAN SYMBOL DESCRIPTION 








BelA2 8,, see figure A-3 
DBETDX 18 


between adjacent streamlines 
dx 


VA2 eee axial velocity 

W2 W., see figure A-3 

CL Axial distance between stations 
WR2 Radial component of velocity 
WU2 Wu2 , see figure A-3 

VU2 Vu2, see figure A-3 

AMR Relative Mach No, at rotor exit 
eZ T 2 

ZS i 

EZ P> 

PRAT 2 [Total-static pressure ratio] * 
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INPUT DATA AND APPROXIMATIONS 


CALCULATE OUTLET ANGLES 
COMPUTE M 


SOLVE EQUATION OF MOTION AND COMPUTE 
STATOR EXIT PROPERTIES 


COMPUTE STATOR LOSS COEFFICIENTS 


ADJUST M, 


COMPUTE ROTOR INLET CONDITIONS 








SOLVE EQUATION OF MOTION AND 
COMPUTE ROTOR EXIT PROPERTIES 








ADJUST M, 





RE-CYCLE ALL OF THE ABOVE 
ACCOUNTING FOR STREAMLINE CURVATURE 


PRESS. 
RATIO 
OK? 


YES 
PRINT RESULTS 


END 


FIGURE A-l: PROGRAM FLOWCHART 
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z CAMBER LINE LENGTH 
=z CHORD 

2 THROAT DIAMETER 

z CURVATURE RADIUS 

z SPACING 

zs MAXIMUM THICKNESS 


TE = TRAILING EDGE THICKNESS PROJECTED IN 
PERIPHERAL DIRECTION 


TN= TRAILING EDGE THICKNESS, NORMALTO FLOW 
- DIRECTION 


FIGURE A-2: BLADE NOMENCLATURE 











FIGURE A-3: VELOCITY DIAGRAM NOMENCLATURE 
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MAIN PROGRAM 


TRAUL 
Sins Pa: XPO 


Y 
ALOS1 


TRAUZ 


= csi 


CSIR XPO 


ALOS2 


TRAUZ 


CSIW 


CSIM Cry 
Gok XPO 


FIGURE A-4: INTERCONNECTION OF THE SUBROUTINES IN THE 
TRAUPEL METHOD 


ALEAK 
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Abb. 8.4.2 Grandwert £5, des Profilverlustes fir Beschleunigungsgitter und Machzahlkorrektur z,, 


FIGURE A-5: INITIAL PROFILE LOSS COEFFICIENT AND MACH NUMBER 
| CORRECTION FROM TRAUPEL 
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| a 
“a anna yo 
“CCT Z 
Pim) yay 
Ga ELE: a 
fa eee 

Q2 , : 05 

7 as 08 Q7 g6 Qs 


Abb. 8.4.4. Korrekturfaktor 73 und Mischverlust ¢ 
infolge endlicher Austrittskantendicke oder 
Ablésung 





Abb. 8.4.5 Facherveriust or 


FIGURE A-6: T.E.THICKNESS CORRECTION FACTOR, MIXING LOSS 
COEFFICIENT AND FAN LOSS COEFFICIENT FROM 
TRAUPEL 
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Untere Werte nur bet Profilen, 
aie weng emptinalich 
auf Zustromwinkelangerung 










Es | 

CCECe 

SRE VE 
SSS 





Sfl—— 


Abb. 8.4.7 Randverlust in Turbinenschaufelungen 


ya a ee ee eee - — - 


FIGURE A-7: "REMAINING" LOSS COEFFICIENT FROM TRAUPEL 
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CALCULATE CRITICAL PRESSURE 
RATIO AND FLOW FUNCTION 


DU FOR EACH S.L. 


<a 


XEe | - (ree ra 





} *FLOW FUNCTION d = CRITICAL FLOW FUNCTION 





FIGURE A-8: 


CALCULATE FLOW AREA 


STATOR ROTOR 
ROTOR 


OR 
STA ADD TIP CLEARANCE 
NO 


COMPUTE FLOWRATE ( M ) 


“ ae 


RAISE VA(3) 


RETURN 


SUBROUTINE FLOWR FLOWCHART 
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CALCULATE ROTOR OUTLET ANGLES 
FOR STREAMLINES 2 AND 4 


CALCULATE = FOR EACH STREAMLINE 


DO FOR EACH S.L. 


SOLVE THE EQUATION OF MOTION 
AT THE ROTOR EXIT (STATION 2) 








NO f 
HAS 
SOLUTION 
CONV. 
r 
YES 
CALCULATE ROTOR EXIT CONDITIONS 
RETURN 


FIGURE A-9: SUBROUTINE ROTO2 FLOWCHART 


a 








0 <0 40 60 €0 90° 100 120 140 160 
4a—~ | ek, 


Abb. 8.4.11 Faktor Ky, fiir Spaltverlustberechnung 


—-_ 


FIGURE A-10: TIP LEAKAGE LOSS COEFFICIENT PLOT FROM TRAUPEL 
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SHLVNIGYOOO SANITINVAYULS 


Piles Veale 
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APPENDIX: B 


DERIVATION OF EQUATIONS USED IN THE PROGRAM 
B-1. EQUATION OF MOTION FOR RELATIVE FLOW: 
The equation of motion for relative flow Ref. [5] is 
—_—> is —— 
Va» WX(VXWraW)+ TVS 
(ee) 
Using cylindrical coordinates, the terms of EQN (B-1) may be 


expressed as follows: 





_ to dHe ,~. dHe .- OH 
VHe - FF @e Foe ae ar or (B-2) 
Le C? tr 
r a 
VXW = in jo 2 
I6 OZ or 
r Nae We Wr 


(9g [oe oWa} be | o(rWa) _ 2Wr 
az or ia or ae 


Lr od Wa _ 2(r Wa) 
r aS ae (B-3) 





a5 





wx (VXW) + 


Co 

















a 7 

— We i 

: OO \N | 3Cr Wi) 7 Wr 1 
ia pia (28 oS -% | +r so 
, sae I (ao _ 2(rw) } - 
[oe (BE = a (2-2 
| 1 f 9Cr Wu) 2\Np ) (20 We | 
(B- 4) 


W xX 200 = ie War + bo We + be Wr [x tg2 


= Ur (2WWa)- bo (QW Wr } 


(B-5) 
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z - 1 Qs es eS 
Tvs = T( r- 338 a C3 2 ee zs 

(B-6) 
Combining equations (B-f) through (B-6) the terms in (B-2) 


can be written as: 











| QHe — Wa] AWe | alrww) | _ 
r ia ae ee 
Wr 9 (rWu) OWr ia S 
r or - ee 6 
(B-7) 
dHe = WW, dW Wa _ Ww o Ve - 
ot a3: its r ae 


(B-8) 


epee ee 
a a aie 
— 


JHR War 3(rWar) aWr 
or a6 = 


o Wr aWe ais 
ry | cee ot Nee <a 
~ ot or | aes We #7 oF 


(B-9) 
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Since the flow has been assumed to be axisymmetric, all deri- 
vatives with respect to 9 are zero. Thus, Equations (B-7), 


(B-8) and (B-9) reduce to, respectively: 


Wa 9 Ur Wun) Wr > (r Wu) — 2W Wr 











- Fr oF  °&4'F or (B-10) 
dHe - Wy dWr = do Wa Wow 
22 a2 ees r 
3 Cane : OS _ 
= 47 a (een 








dHe _ Wu 3CrWa) _\, We t We We 








an 
or Tr er o> or 
aS 
bs War + 
2 wand ie (B12) 
Equation (B-10) may be written as 
> Cr Wa) a! Wr Oo Ur Wu) =) 2900 ¢ Wi 
gz 7 Wa or We ~~ (B-13) 
Substituting into equation (B-1l), 
dHe Wa Wr alrWe) , wy, Wr 
o% ro Wa or oz 
Or Wee : Q 2 (B-14) 
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Multiplying equation (B-9) by Wr and Equation (B-14) by Wa 
meciits In 


w, ote = WuWr 30rWo) _ ww, = 








or r or . 
We. Wir a aT 20 Vw Wy t Wr T a (B-15) 
and 
gHe ~) WaWr alrWw) Wr 
Wea A — " -, — + ee = 
We Wy; —2We 2 WW Wr + Va T ae 
o- a2 
(B-16) 
Adding these two equations yields 
NV dHe Wy S, ore = = T|W VES oo 
or dz (B-17) 


Since the flow has been assumed to be adiabatic, the total 


Belative enthalpy, H 1s constant along a streamline. Thus, 


Re 


aHR dHR 
Vin oe dz bid oF (B-18) 


and re-arranging, 


aoHe a. Wr d He 
92 We or 
(B-19) 
From equation (B-19), eq. (B-17) can be written as 
aS Wr aS 
ee (B-20 
at Wee ar ; 
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puboereucing Eq.s {B=19) and (B-20) into equation (B-15) gives 








_ We dhe | We Wr 3CrWue) yy 2Wr 
We of r Wo. or ee 
7 Wo. Wa Wi; Ww 2s 
We = -2Qw —eOr LL r 
: or We Wa q or 
(B-21) 
foltiplying Equation (B-21) by 7-2 yields 
R 
HR = Wa 9 (rWu) yy Wr _ — 
or r or ~ 32 
ol WW 
Wa — rQwWeteT 2S_ ee) 


or 


This expression is identical to equation (B-21) and is the 
equation which must be solved. It must be put into a form 
which can be solved by the computer. Re-writing equation 


e-22)) siven that 





o Nev bells na 2 Cw,”) 
—— so ar 
BiG.) _ oWr 2 iw 3 Cr War) 
yields - or! aAWa a> rt - or T 
UO \W ae dHe as — O 
4 eS +t ade 


The relative enthalpy can be written 


Pee We 
age J X46) z 2g eV L4eU 
(B-24) 
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iiemcamrvalent enthalpy, defined in ref. [1] is 





> q 
He = i + wt + ae 
qe ac (3525) 


emcee, the nelative enthalpy can be written as 


Zz 
= Us 
He = He - > 


(Be 26) 


Also, the turbine outlet static temperature can be written 


aS 


He Wa" 


ay 


P ACp 





(B-27) 


substituting Eq. (B-26) and Eq. (B-27) into Eq. (B-23) and 


Beery ing Eq. (B-21) to the rotor exit, gives 


_9 CW") = 2Wa a Wr. i, 5) Warr 3. Cr Wire) 








ae 
oh z a? r ole 
es 
4 w Wu, - — He- = +2 ee ee 
Cp 
oe 
Wa dSa 
= | or, —° 
= (B-28) 
erven the relationships: z 
Tan? y= 
PN = We 





| 
ol. es eee 


Equation (B-28) can be written as 


_dCWa; ) — 2\Vo AWrs _ War 


oz = Be Cp Cos fee. 5. 
Wu2 2 (Wu) 
= — +4 Win, - o> oes He. 


ere.) ot cae 6 
o Vp. . Cp tHe — Wa 





a 
a 
and substituting 2 Us") = aD 


gives 


o a. | 


Asa 
oT, 


(B-29) 


into equation (B-29) 
ole 


men .3(3rororr eS ons 
ae. be 


2 Wawa 2 Cre Wa.) 


HE 
- xe + YW Ww,-2k 





Sa 
W. 
QU" eZ Late a a 


Multiplying Eq. 


a 
m 


less form of Equation (B-29): 


~ 


= 50) 


G 
(B-30) by a) results in the dimension- 
W 


a 





Yam _ 2 CWa,) 6 5 War Wor 


- 
We, i Ge We, War “™ 32 
[2 Wate War 
Wee Yom aise 2) Wawa Waa Tom 3 Tom Warm Wa. 
a s nT eee 
Warn Cpe a2 = A Wa, Var V2 d (re / Tem) 


py lam Wir Was 5 tam dre, Witte Fem _ 


- Pe 2 + 
We, an Waar MoU ae on VMN 
~ al 
Tam Ate _ Wares Was Sr 
2 a ——- ae 
Co ec Was. vaee. ae 
(Bao) 
Introducing the non-dimensional quantities 
a: Wa. (B-32) 
Warm 
Co mn 
YX = = a 
: r 
i (B-33) 
¥ S (B-34) 








Pomation (8-31) iS written as 








30), Y' - awe  _y _3s* 
OX We. a2 cos") OY 
TanS_ 9 (XY Tan6) Um YTANQ 

eT x aX _ Was en 

as oHeE ie Uy aHe 
a —- rer/? ll Oe —— 
Wo, oY es Wa'm 

2 a he _ 

Y" Tan pi n=. we, 


—, 


Mrs fourth term of Eq. (B-35) is 


2Y os d ( <etows) r 2Y a X % 





a te) 


d Tan & T oY FF dX 
meee *TOND +t Tene Sy 


= 2Y° TAN ¢ TanB 22 OY 
8 —— + 2Y Taw 6 ae + 


a 


Y 2 
a x Tan 6 
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also, aaa R _ B55 
3X Cos" aX 
“LY Tat @ 2 = Tan @ oe 


equation (B-35) can be written 

SE 

9) (14 Tan’ )- 3 bie r, 2WVr _ y"_ as" + 
o ¥ Wa 2 cos? OY 
YUmY TANB _ 


2 Tan® ao ee 
x Wee wr. 


A7 cos*B OY 


Therefore, 














2 aHe Um U $ 
we oe? ae | Wa THO Se 
Bonn Wa mn Wie oy 

(ba56) 
cos” G 

Multiplying Eq. (B-36) by ——a ) and observing that 

( + TAN® 8 = ! ny 

s*A 


me 3Cy*) 2 2 er 
TO OSE WOR” Cok Ok 


ye 
4 Um Sin B cos B ‘ a Um Us co¥’B 
ah wal 


7 oh Zz 2 
& IANB + —— sin Bt 
o X % P Wa, ¥ Wain Y 


/ 2088 QHeE 2HeE cosB 2 2<* 
(B- 37) 


5 : 
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iom@eecount for streamline curvature the following terms are 
introduced: 


me 


2 AR \* 
' + (2) 


Cos” ) = 


(Basen 
where A, the angle between the axial and radial components of 
velocity at a point, is approximated as the average value be- 
tween two stations. 


Aso, 





kK &R | aWr 
a Wa oh 
(B- 39) 
where 6r is the streamline shift throught the rotor defined 


as Sr = Taotor — ‘ROTOR 


CUleoe t Nine ws (B-40) 
Substituting Eqs. (B-38) and (B-39) into (B-37) yields 


Miley) = 2,\_ SR L7+ (ary | ds? 
Tay = cos’ 8 (2K Cm <a |) = ( a er 


E d8 2 2 4 Um SIN B 
OMAN) se SIN _ 1Um SING cos B ee 
arr p Wa Y 


Xx 
OW pr 


2 Um Va cos? 2cs 8 due _| ane cos'B 22] ols* 
Wee 2 a Wary? dy we. v2 ~SInB AY 

On any Q im a, 
(Beane) 
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To obtain a dimensionless equation, the term 
eS = AGe J 


is introduced into Eq. (B-41) giving 


d(In¥*) _ aal_y alsrirm) (Lae) ds* 
ate). — conta (x Heep). [e ene | de 


ata 








2Tan - - af siN7@ - 4Um Sin bcos G = 2Vy Uo cvs" 
x * Wan Y Wa? Y> 
Cy cos ® 8 dHe eee cos 6 _ sin2Z as* 
O a 
War YY? dX . | Wa, Y* a 
(B- 42) 


Equation (B-42) is the form of equation of motion used in the 


computer program. 


B-2. EQUATION OF MOTION FOR ABSOLUTE FLOW 


The equation of motion for absolute flow 


@baaS) 
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Differs from the equation of motion for relative flow 


VHezwx(vxw +aw )+TVS 
- (B-44) 
only by the term Wx2w which is the Coroiolis acceleration. 
To obtain the programmed form of the eauation of motion for 
the stator, the previous derivation is followed, but with 


Ne= 0, HH, becomes H, W becomes V, and 8 becomes a. 


foo eit AREA RESTRICTION FACTOR Z 

The condition at the outlet of a blade row with boundary 
layers on both sides of the flow channel is shown in Fig. B-l. 
The flow 1s considered to be turbulent within the boundary 
layer while, outside the layer, the velocity of the flow is 
the theoretical velocity. Assuming a power-law velocity pro- 


file, the velocity may be written, 


Been e . 
rH 8 


Se ar | (B-45) 


The mass flow rate exiting the blade row can be expressed as 


Z 5 s§ 
2 Wenicos & ee oo Oe ee 
- . aig ae » S COSA) Cos XX 4 


2) awedy 
where O,. and en represent the ideal conditions for an 


isentropic expansion through the blade row to the discharge 
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pressure Pd, which is assumed to be constant across the blade 
spacing. The discharge angle of the flow leaving the blade 


row is closely approximated by the expression [Ref. 1] 





= OW 
Kg = Cos + 
a 
cose y (3-47) 
Inserting Eq. (B- 47) into (B- 46) and reducing yields 
& 2 ML 
Mm =P ry Vey Oo a (| , (tH an 
(B-48 ) 
Assuming a perfect gas 
fp TTH Tro — (Tre - Tra ) 
i nase eee 
(rH i Tro - ae ‘ra 
(B-49) 


Defining = 


Pa 
ee =l- ( ye (B-50) 


Equation (B- 49) can be written 








e |- Xe 
—————— Zz 
Me | 
| Hr }- Xe (=) (B-51) 
PmMestecucame Eq. (8-51) into (B- 45) yields 
$ 4K" 
t= Civ 2 |l-2 (I(x) | 2 dy) 
0 


(ebo2) 
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Using the displacement thickness given by 


Vr 


g "26. (140) (78a ol h 


(B-53) 
the mass flow ratecan be written as 
* 
re 
IM = Pra VTH OL 
(B-54) 


The loss coefficient, expressed in terms of average kinetic 


energy lost is 


aE E 
‘ = <a o = | — ; a 
mM ( Ym) Wr TH 
> od 


(oso) 
Poreres =: 15 the actual kinetic energy of the flow, given by 
y, - aw 
TH 
seminars) (pe 
O 


(B- 56) 
Bemipstitucimea Eq. (B-51) into (B-56) gives 
a 3m 
F-@, Ve ££ (1-U-xe)| tay 
a (1- Xen”) 
(B-57) 


dak:0 





The energy thickness is written as 


| 3m 
¥* a e 
§ = § ie (1-Xe ) | ( \ 2m dh 
o (i Xe”) 
(B-58) 
The loss coefficient can therefore be written as 
a 
é | | — Z OU 
= $* 
= 9 B- 
I- 2 (B- 59) 
The area restriction factor 2Z, is the fraction of the flow 
area through which the uniform theoretical velocity would 


produce the actual flow rate, thus 


aS” 
_ OL 





(B-60) 
Defining the energy parameter (a form factor) as 
we KL 
ot § 
= oo 
3” 
Esa on) 
using Equations (B-59) and (B-61), Eq. (B-60) becomes 
yy eeK 
aa KY 
Ho" =-1 + Pp 
(B-62) 


where €p is the profile loss coefficient. 


ee 





Sere | LHE ENERGY PARAMETER, ea 
In Equations (B-53) and (B-58) the denominator of the 


integrand is expanded using the binomial theorem, so that 


2n\-7! of In os 
(1- eh 2 | +Xeh + Ye A 7 Xe he. 93 


The integral of Equation (8-58) is now written as 


, a“ | mn m Pn 
( a ah ( 1’ ae » Kern : 











_ 2m 
o I- Xe 2 
3 Im yo ilm 
Len teh | AA (B64) 
which, on integration becomes 
3m 2 
5 
| ee al n = + ee ee tee 
Im Xe/ ee 3mt+l Smt “~mH = Amt 
. / (B- 65) 
Him! ,?! 
Therefore, Equation (B-58) becomes. 

KK x Ve Xe 
ee te a se, 
S 3m +! Smt “m4 Gm +! 

Ye ( x.) (B- 66) 
Il mt I Kg 








which can be written as 








#4 aL 
& Xe-l 3m+l Smt ]m +1 


2 4 
Xe, Xe! 
Tm+l Il mt 
(B-67) 


In a similar manner, 





a ¥ 2 3 
a (Xe -1) ———- + pe eS 
§ Xe -l Mtl sm4 Smt “)m tl 
L 
XE 
5 Sim +1 


(B-68) 
Substituting Eq. (B-67) and Eq. (B-68) into Eq. (B-61), the 


equation for H*** used in the computer program is obtained: 








a 3 1 
| | ie x XE NEL 
= i fa + aay 
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Eee NU 
PROGRAM SEGMENTATION ON THE HP-1000 


Segmentation allows large programs to be run on the 
neewo00. The program is divided by the programmer into a 
main program and several segments, which are stored on the 
disc. Each segment and the main program are then compiled 
and loaded. When the program is executed, the main program 
and its segments are called into memory individually, and 
Mitiyveasstney are needed for execution. In this manner, a 
program can run in a partition which is smaller than that 
program's total size. 

When the main program has performed all executable state- 
ments, the first segment is called into memory by an EXEC 
call. The system then loads that segment from the disc into 
a memory block following the end of the main program. The 
meeoecess 1S illustrated in Figure C-1. Note; the main program 
plus the largest segment may not together exceed 29 k. Once 
fPesocoment 1S in memory 1t can call another segment. 

When executing, any segment can call any subroutine which 
is attached to the main program. It was this feature which 
allowed the present program to be run. All subroutines were 
placed within the main program. In fact, the main program 
consisted of nineteen subroutines and functions. A segment 


may not return to the main program. Communication of data 
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between the main program and the segments is accomplished 
through a common block. 

The four segments of the present program are ''MAIN", 
"SHORT", "PART 2' and "PART 3'. The manner in which control 
is passed from the main program to the first segment and from 
the first segment to the second is as follows: 


BLOCK DATA 


END 

PROGRAM THESS 

DIMENSION INAM (3) 

DATA INAM /2HSH, 2HOR, 2HT / 


CALL EXEC (8, INAM) 

END 

PROGRAM SHORT (5) 

DIMENSION INAM (3) 

DIMENSION NAME (3) 

DATA INAM /2HSH, 2HOR, 2HT / 


DATA NAME /2HPA, 2HRT, 2H2 / 


CALL EXEC (8, NAME) 
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END 
PROGRAM PART 2 (5) 
DIMENSION NAME (3) 
DIMENSION NAMR (3) 
DATA NAME /2HPA, 2HRT, 2H2 / 
DATA NAMR /2HPA, 2HRT, 2H3 / 


CALL EXEC (8, NAMR) 

END 

PROGRAM PART 3 (5) 

DIMENSION NAME (3) 

DIMENSION NAMR (3) 

DATA NAME /2HPA, 2HRT, 2H2 / 
DATA NAMR /2HPA, 2HRT, 2H3 / 


END 
The "(5)" after the program name indicates that it is a 
program segment. Note the manner in which the program name 


is put into a data statement using the Hollerith notation. 


a7 











SEGMENT l SEGMENT 
SEGMENT OVERLAY 


AREA 


+e 


MAIN 
PROGRAM PROGRAM 
AREA 


PARTITION 
AREA 


SEGMENT 2 





DISC MEMORY 


SloteM 
TABLES 


LOGICAL MEMORY 


FIGURE C-1: PROGRAM SEGMENTATION-ILLUSTRATION OF THE MAIN 
PROGRAM CALLING A SEGMENT INTO LOGICAL MEMORY 
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APPENDIX: D 
RUNNING THE COMPUTER PROGRAM 


If the reader is unfamiliar with the HP-1000 Computer 
System, references [11] and [12] should be consulted before 


attempting to run the program. 


D-i. DATA INPUT 

Using the editor, input the following data into segment 
SSoHORT'. 

1. Turbine operating conditions: referring to Table 
A-III, type in appropriate data in lines 66 through 69 and 
74 through 78. 

2. Special input data/program control parameters: 
referring to Tables A-IV and A-V, type in appropriate data 
in lines 83 through 98. 

3. Turbine geometry: referring to Tables A-1 and 
A-II, type in data for stator and rotor in lines 103 through 


1S 6 . 


D2 COMPILING THE PROGRAM 
1. To compile the main program type: 
eee INe ,.MAIN::25,-,- 
2. Compile the first segment: 


:RU,FTN4,SHORT::25,-,- 


ETS 








Secempmte Che second segment: 
eReeNa ge ART 2::25,-,- 
4. Compile the final segment: 


TRUE PARI S:°325,-,- 


LOADING THE PROGRAM 
iyasie 
RU, LOADR 
Ham return key 
Will display 
LOADR: 
Type 
OP,LB 
Will display 
LOADR: 
iyve 
SRE, 3MAIN: :25 
Will display 
LOADR: 
Type 
gee oomMORT::25 
Will display 
LOADR: 
Lype 
MEMO ARI 2:2 25 


Will display 





LOADR: 

ype 
PREG SPARE Os. 25 

Will display 
LOADR: 

Type 
> END 

Arter the end statement, the loader will display that the 


program is ready for execution. 


D-4 RUNNING THE PROGRAM 

Type 

RUN, THESS 

The program will be executed and no further action by 
the operator is required. The computed pressure ratio of 
each iteration of the outer loop of the program is displayed 
on the screen as it is calculated. The operator therefore 
has some idea where in the iteration process the computer 


pRooram 1S executing. 








AE BeEN DL XS EC 


DISCREPANCIES IN MACCHI'S PROGRAM 

1. Main program, lines 21 and 22; the value of ICL has 
not yet been read. 

2. Main program, lines 163-166; the Traupel method of 
calculating gas outlet angles does not take the Mach number 
into consideration. However, in lines 163-164, the program 
is attempting to draw a parabola through points which repre- 
Selle Outlet angie as a function of Mach number. 

3. Main program line 281; the calling of subroutine 
SLINE is questionable. Parameters are transferred to that 
subroutine, but many of them have not yet been defined 
(HE, DHEDX, WPER2, DSDX1). These undefined variables will 
be set equal to zero by the IBM 360 and 370 computers. Thus, 
in line 10 of subroutine SLINE, the value of DWDX will be 
zero and in line 17, division by zero will occur and the 
execution of the program should cease. 

4. Subroutine ROTORI lines 22 and 26; the stator radii 
are used in the calculation whereas the rotor radii should 
be used. 

S. Subroutine ASOSI, line 107; the correct Fortran code 
1s 


ZemeePS(l) = .5 * ZETAS(T) 


eZ 





6. Subroutine ALOS2, line 121; the correct Fortran code 
1s 
GEGAPR(T) & .5 “S2ZETAPR(T) 
Peeoupbroucine ALOSZ, lines 123-126; the stator radii are 
used in the calculation whereas the rotor radii should be used. 
8. Subroutine ANGAIN, line 14; the correct Fortran code 
1s 
AOp=PATAN(T. -XCL/H*CH*C0S (ANGI) /COS (ANGZ) * 
TAN (ANG2) + XCL/H*CL*COS (ANG1) /COS (ANG2) * 
TAN (ANG1) 
Note: Since reference [2] was published, Professor Macchi's 
program has been further developed by Professor Macchi under 
private sponsorship [Ref. 13]. The new code however, is not 


generally available. 
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APPENDIX F 


COMPUTER OUEPUT 
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